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Two basic physical models, a two-level system and a harmonic oscillator, are realized on the mesoscopic 
scale as coupled qubit and resonator. The realistic system includes moreover the electronics for controlling the 
distance between the qubit energy levels and their populations and to read out the resonator's state, as well as the 
unavoidable dissipative environment. Such rich system is interesting both for the study of fundamental quantum 
phenomena on the mesoscopic scale and as a promising system for future electronic devices. 

We present recent results for the driven superconducting qubit - resonator system, where the resonator can 
be realized as an LC circuit or a nanomechanical resonator. Most of the results can be described by the semi- 
classical theory, where a qubit is treated as a quantum two-level system coupled to the classical driving field 
and the classical resonator. Application of this theory allows to describe many phenomena for the single and 
two coupled superconducting qubits, among which are the following: the equilibrium-state and weak-driving 
spectroscopy, Sisyphus damping and amplification, Landau-Zener-Stiickelberg interferometry, the multiphoton 
transitions of both direct and ladder-type character, and creation of the inverse population for lasing. 



Contents 



1. INTRODUCTION 



1. Introduction 



[11 



2. Semiclassical theory of the qubit-resonator system |3 

Krylov-Bogolyubov formalism for qubit-resonator 
system 

Inductive coupling with LCR resonator. Parametric 
inductance 

Capacitive coupling with nanomechanical resonator 
Parametiic capacitance 



3. Dynamical behavior of a two-level system 

4. Excitation of a superconducting qubit 

Inductance of superconducting qubits 

Equilibrium-state measurement 

Resonant transitions in the charge qubit 

One- and multiphoton transitions in the flux qubit 

Interferometry with nanoresonator 

5. Multi-qubit systems 

Equations for a system of coupled qubits 
Weak-driving spectroscopy 
Direct and ladder-type multiphoton transitions 
Lasing in the two-qubit system 

Conclusions 

Acknowledgments 

References 







12 



12 
12 
13 
H 

M 



18 



A quantum system, subjected to external driving, can expe- 
rience resonant transitions between its energy levels. Conser- 
vation of total energy assumes absoiption or emission of sev- 
eral photons of the driving field. Such multiphoton processes 
play an important role in atomic and molecular systems inter- 
acting with electromagnetic field [1]. For example, the multi- 
photon resonant spectroscopy is one of the methods to probe 
the structure of atoms and molecules [2]. This technique has 
the advantage of observing highly excited states by using rel- 
atively low frequencies. The concept of another application, 
the multiphoton excitation microscopy, is based on the mul- 
tiphoton excitation of the fluorescent dyes molecules 13-5]. 
This technique allows imaging biochemical objects with high 
spatial resolution. 

Recent development of fabrication and measurement tech- 
niques enabled a study of the wide spectrum of quantum phe- 
nomena in superconducting structures. During the past years 
it has been clearly shown that specially designed macroscopic 
superconducting circuits, which include Josephson junctions, 
behave quantum mechanically similar to a quantum particle 
in a potential well. Under certain conditions, these objects 
demonstrate the coherent superposition between their macro- 
scopically distinct quantum states. It is important to note that 
this is a pure quantum effect which has no classical analogue 
and can be used for a number of intrigued applications. If the 
circuit's dynamics can be described in the frame of the two- 
level approximation, such two-level quantum system is called 
a qubit. The advance in the study of different phenomena in 
superconducting qubits can be found in the reviews 

In general, superconducting Josephson circuits can be de- 
scribed as multilevel quantum systems. By analogy, such sys- 
tems are called aitificial atoms, while coupled qubits systems 
behave as artificial molecules. An interesting problem is how 
phenomena, known from atomic physics, will appear for these 
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artificial atoms and molecules. Note that the following fea- 
tures differ these mesoscopic-size quantum systems from their 
microscopic counterparts: a high level of controllability by 
electronic means, coupling to the macroscopic-size read-out 
devices, and unavoidable dissipative environment. 

For characterization and controlling the states of supercon- 
ducting qubits the one-photon spectroscopy was done by us- 
ing relatively weak driving Lli-JJlJ . Matching of the ground 
and higher states with the one-photon energy was exploited 
to probe the upper levels of the Josephson-junction circuits 



1 19l - l26ll . With increasing driving power, the multiphoton ex- 
citations were used to study the features of the artificial atoms 
both for the two-level dynamics II27H3III and when the upper 



levels were involved Il32l - l36ll . For strong driving, the width 
of the resonance lines periodically tends to zero, which can 
be described as the destructive Landau-Zener-Stiickelberg in- 
terference [37]. Respective interferograms displayed double- 
periodical dependence of the upper-level occupation probabil- 
ity on the energy bias and the driving amplitude [i38M43il . 

Two and more coupled qubits can be treated as artificial 
molecules. Being excited by a resonant microwave field, 
they display one-photon transitions 1 441 45 ill . Alternatively, 
at smaller frequencies, the two-qubit systems can experience 
multiphoton transitions 1152145 511 . 

In this article we review the observations of the multipho- 
ton transitions in single and coupled superconducting qubits 
probed by a classical resonator, and also we present the re- 
spective theory. Having the purpose of presenting and de- 
scribing specific results for the multiphoton transitions, our 
consideration is limited to the Josephson-junction qubits. We 
note however that similar phenomena can be studied in dif- 
ferent quantum objects, which can be described as two- or 
multi-level systems, such as quantum wires and dots I5^60ll. 
nitrogen vacancy centers in diamond llil , 62], ultracold atoms 
||63-65|, nanomechanical and optomechanical setups 166 -68.1 . 
electronic spin systems, two-dimensional electron gas, and 
graphene I69l-7lll. 

The paper is organized as follows. In Sec. 2 we use the 
method of an asymptotic expansion for the qubit-resonator 
system in order to obtain the resonator characteristics. This 
formalism allows us to separate the dynamics of the relatively 
slow resonator and fast qubit. Then, in Sec. 3, we consider the 
multiphoton dynamics of an isolated two-level system. Later 
the formulas of those two sections will be applied for the de- 
scription of the experimentally observed multiphoton excita- 
tions in single qubits (Sec. 4) and in coupled qubits systems 
(Sec. 5). 



2. SEMICLASSICAL THEORY OF THE 
QUBIT-RESONATOR SYSTEM 

For characterization of a quantum system different tech- 
niques can be applied. One of the possible solutions is to use 
the so-called parametric transducer [72 1. A key element in 
any parametric transducer is an optical or a radio-frequency 



auto-oscillator. A transducer, coupled to the quantum system 
of interest, is constructed so, that quantum system dynamics 
causes a change of the phase or/and the amplitude of its oscil- 
lations. A phase (amplitude) shift provides information about 
the dynamics of a quantum system. In particular, for prob- 
ing the qubit's state, several types of oscillators have been 
akeady used: an LC tank circuit [73, 74], a nanomechani- 
cal one [I43, 75 1, and a transmission line resonator (76, 771. 
If the resonator quantization energy hiOp is smaller than the 
thermal excitation energy k-^T, the resonator can be consid- 
ered as a classical oscillator Then the qubit-resonator system 
can be treated semiclassically: here a qubits quantum system 
is driven by a classical field and probed by a classical oscil- 
lator It is important to note that the similar approach is well 
known in quantum optics - many phenomena in the atom-light 
system can be described by making use of this semiclassical 
model litl. 



In this work we present the semiclassical description of 
some observed effects for the resonator-qubits systems. We 
will not consider here the situation of coupling the qubits sys- 
tems to a high-frequency resonator, which can be realized as 
a transmission line resonator The quantum properties of this 
qubit-resonator system are not described by the semiclassical 
model. For recent works in this field see e.g. [78-82] and ref- 
erences therein and also Refs. [js^-El. where the multiphoton 
excitations were used to drive transitions between the multi- 
ple energy levels of the qubit-resonator system in the strong 
coupling regime. 



Another note here should be made about the term "mul- 
tiphoton processes". In the context of the semiclassical ap- 
proach, it relates to the energy of several photons which is ab- 
sorbed or emitted by the quantum system. In the broader sense 
the term "multiphoton" can relate to other processes employ- 
ing the quantum nature of the electromagnetic field, see Ref . 
[86] for a review of the non-classical phenomena in entangled 
multi-photon systems. 



This section is devoted to the properties of the qubit- 
resonator system. It will be shown that in the frame of the 
semiclassical approach the influence of the qubit on the res- 
onator can be described by the "renormalization" of the os- 
cillator constants. For instance for a mechanical resonator it 
can be quantified by introducing the equivalent qubit's-state- 
dependent elasticity coefficient and damping factor In the 
case of inductive/capacitive coupling, the qubit's impact on 
the resonator can be described by introducing the qubit's- 
state-dependent effective inductance/capacitance, while the 
losses can be described by the effective resistance. For con- 
creteness, we will consider two realistic systems: the flux 
qubit inductively coupled to the tank circuit iIstIi and the 
charge qubit capacitively coupled to the nanomechanical res- 
onator [881I . 
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FIG. 1 : (Color online) Qubit (quantum two-level system) coupled to a classical resonator, (a) Schematic diagram of the model qubit- 
resonator system. The qubit is represented by the two-level system with the two states, |— ) and |+) , and with the energy difference /S.E. The 
resonator is demonstrated as the spring oscillator with the elasticity coefficient fco. As described in the main text, influence of the qubit on the 
resonator can be described by introducing the effective elasticity coefficient fceff , which includes the qubit's-state-dependent (or, parametric, 
for brevity) elasticity coefficient fcq. (b) The flux qubit coupled via the mutual inductance M to the LC resonator. This can be described by 
introducing effective qubit's-state-dependent inductance L^a, which includes the parametric inductance Lq in parallel to the tank's inductance 
Lq. (c) The impact of the charge qubit on the nanomechanical resonator's state can be described by introducing the effective qubit's-state- 
dependent capacitance Ccft, which includes the parametric capacitance Cq in parallel to the resonator's capacitance Co. 



Krylov-Bogolyubov formalism for qubit-resonator system 

First let us consider the mechanical resonator as a spring 
with the elasticity fco, the damping factor Ao (which is as- 
sumed to be small), and loaded with mass to, as shown in 
Fig. [Ila)- The oscillator has eigenfrequency wq — ^Jk^jm 
and the quality factor Qo = "^'^o/'^o- Its state is influ- 
enced by the qubit through the force eFq and is driven by 



the probe periodical force eFp sin Wpt. Here the small param- 
eter e is introduced explicitly to emphasize the small qubit- 
resonator coupling as well as the amplitude of the external 
harmonic force eFp, which enables us to make use of the 
asymptotic expansion method. The external nonlinear force 
is assumed to depend on the variable x and its derivative only, 
Fq = F^{x,dx/dt). 



The displacement x is the solution of the motion equation 



(fx 



dx 



TO-^ + Ao— + ki^x = ei^q ( a;, — ) + eFp sinwpt. (1) 



dt 



dx 



dt 



The oscillations in the nonlinear system described by Eq. ^ 
can be reduced to oscillations in an equivalent linear sys- 
tem by making use of the Krylov-Bogolyubov technique of 
asymptotic expansion ll89ll . Specifically, in the first-order ap- 
proximation with respect to a qubit-resonator coupling pa- 
rameter and close to the principal resonance, Wp sa ujq, 
the equivalent linear system is characterized by the effective 
amplitude-dependent elasticity coefficient kcs{v) and the ef- 
fective damping factor \cs{v) (see chapter 7 in Ref. ll89ll ): 



'df^ 



dx 

Acff(w)-^ + kes(v)x 



e-Fp sinojpt, (2) 



V cos(ajpi + 5), 



(3) 



2n 



fco / Fc^^iv, ip) cos = fco + fcq, (4) 



= Ao + 



TTVUJQ 



i^q(w,-0) simpd^p = Ao + Aq(5) 



where Fq{v,^) = Fq 



dx ^ 
dt . 



Fq {v COS ijj, —ujpV sin -0) 



Note that in Eq. (|2]i both v and S are time-dependent values. 

Here, in equations ^ and (|5]l, we have introduced the para- 
metric elasticity coefficient fcq and damping factor Aq. In this 
context the adjective quantum is sometimes used instead of 
"parametric" to emphasize that it is the qubit-state-dependent, 
i.e. it is defined by the quantum properties of the coupled sys- 
tem. In what follows, by simply changing the notations we 



will see that the parametric elasticity coefficient gives either 
parametric inductance or parametric capacitance, when cou- 
pling is inductive or capacitive respectively, while the para- 
metric damping factor will give us the parametric resistance. 
Note that in equations ^ and (|5]l the parametric terms kq and 
Aq are of the first order in the small parameter of the problem 
e. 

This linearization procedure allows to obtain important in- 
formation even without solving equations of motion. In par- 
ticular, the effective resonance frequency of the linearized sys- 
tem Weft = \/ kcs I m gives the expression for the frequency 
shift 

k 

Alo = Wcff - Wo = -—2— . (6) 
ZmujQ 

For physical interpretations it is important to emphasize 
that the application of the linearization technique resulted in 
the substitution of the nonlinear force by the linear one: 



^ , dx 



dx 

Fq = —kqX — Aq-^- 



(7) 



This latter "parametric" force describes the work done by 
the quantum system over the resonator; the respective energy 
transfer during one period is the following 



2-K/ujp 

W= I Fc 



dx 
'dt' 



dt 



-TTUJnv'^Xn. 



(8) 



This, in dependence on the sign of the parametric damping 
factor Aq, describes periodical extraction or pumping of the 
energy by the quantum system out of or into the resonator 
This is known as the Sisyphus damping and amplification 



The solution of equation ^ in the first approximation in e 
is given by the expression ^ with the amplitude v = v{t) 
and the phase shift S = 5{t) slowly varying in time. For these 
values the asymptotic expansion method gives the following 
system of equations (see chapter 15 in Ref. |89]) 



Acff(w), 
2m 



dv 
'dt 
dS 

— = U}cS{V) - UJp 



■ cos 5, 



■ sin 5. 



(9) 



(10) 



m,v{ui{) + Wp) 

In the regime of stationary oscillations: dv/dt = d6/dt = 0, 
and we obtain equations for the amplitude v and the phase 
shift S, which can be written in the form 



tan 6 = 



LOoXcsiv) ' 
eFp cos 6 
woAoff(u) ' 



(11) 



(12) 



In what follows it will be demonstrated that the phase shift 
S and the amplitude v can be directly observed experimen- 
tally, which gives the information about the quantum system 
through the values of the parametric elasticity coefficient kq 
and damping factor Aq. 




FIG. 2: (Color online) Flux qubit coupled inductively to an LCR 
(tank) circuit. The flux qubit is pierced by the magnetic flux 
induced by the current in the controlling coil and by the current in the 
tank's inductor. The qubit is coupled via the mutual inductance M 
to the tank circuit. The resonant tank circuit consists of the inductor 
Lo, capacitor Co, and resistor Rq; the circuit is biased with an RF 
current /bias- The tank voltage V is the measurable value. 



Inductive coupling with LCR resonator. Parametric inductance 

Now we consider as an illustrative case the system of a 
flux qubit (with geometrical inductance L and average current 
/qb) coupled inductively to the LCR tank circuit, as shown in 
Fig. 121 The approach, presented here, is the development of 
the theory in Refs. ll9ll - l93ll . The quantum system is consid- 
ered to be weakly coupled via a mutual inductance M to the 
classical tank circuit. The circuit consists of the inductor Lq, 
capacitor Co, and the resistor Rq connected, for the specifica- 
tion, in parallel. The tank circuit is biased by the current /bias, 
and the voltage on it V can be measured. 

The flux qubit can be described by the pseudospin Hamil- 
tonian [I93l 



H = (Tr (Tr, 

2 2 
£{t) = eQ + AsiiiLut, 



(13) 
(14) 



where the diagonal term e is the energy bias, the off-diagonal 
term A is the tunneling amplitude between the wells (which 
corresponds to the definite directions of the current in the 
loop) and ax.z are the Pauli matrices. 

To obtain the equation for the tank circuit voltage, we write 
down the system of equations for the current in the three 
branches, namely, through the inductor (/l), the capacitor 
(/c), and the resistor (/r) (in particular, for systems with su- 
perconducting elements see e.g. Ref. 19511 ): 



/bias 

V 



Il+Ic+Ir, 
CoV, In - V/Ro, 



(15) 

(16) 
(17) 



where $0 is the flux through the inductor of the tank circuit. 
This flux is the response of the quantum system to the flux, 
induced in it by the current /l. It follows that the voltage 
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V in the current-biased tank circuit (/bias = /Asinwpt) is 
described by the following nonlinear equation 



Co 



(PV 1 dV 1 



(18) 

The external flux $c is assumed to be proportional to the cou- 
pling parameter = A'P / LLq ^ 1 and to depend on time 
via the voltage V and its time derivative V. Equation (fTsT i for 
the voltage V coincides with the nonlinear equation ([T]i for the 
variable x with obvious change of the notations. 

Thus, the formalism presented in the previous subsection 
is directly applicable for the given problem. Specifically, in 
the first order approximation with respect to the coupling pa- 
rameter k'^ and close to the principal resonance (up ^ uo = 
1 / ^/LqCq), the equivalent linear system is characterized by 
the effective resistance Rcs and inductance L^s as following 



Co 



^cK^ + LJV = lAOJp COStJpt, (19) 



V = V cos{ujpt + 6) , 
1 11 



Rcs{v) Ro Rq{v) ' 

1 11 



ioff(w) Lo Lqiv) 



(20) 
(21) 

(22) 



Here Qo = ojoCoRo is the quality factor of the unloaded tank 
circuit (at $c = 0) and the parametric (qubit's-state depen- 
dent) resistance Rq and inductance Lq are given by the for- 
mulas 



Rq{v) 



Qo 

TTvRn 



$o(w, V') s'mipdip, 




2tt 



TTvLo 



<l'c(w, ip) cosipdip, 



(23) 



(24) 



where $0(^,-0) = ^0(^7^) = <i>o(w cos i/", —wwp sini/j). 
The resonant frequency Uctf becomes amplitude-dependent 
and is shifted by 



which can also be rewritten alternatively in terms of the ef- 
fective quality factor Qcb — <^oCoRcs{v) and effective fre- 
quency shift = [woff (u) - ujp] /uJo- 

Thus, the observable values - the amplitude v and the phase 
shift S - are defined by equations (|26] |. which depend on the 
response of the measurable system, ^e{V,V). As we dis- 
cussed above, strictly speaking, the dynamics of the tank cir- 
cuit has to be considered jointly with the dynamics of the qubit 
(corresponding calculations see e. g. in Js^). However, in 
what follows we consider two illustrative limiting cases, when 
the dynamics of the qubit can be treated separately from the 
dynamics of the tank circuit. For simplification we introduce 
phenomenologically the relaxation time Ti which is caused 
by the coupling to the environment and to the tank as well. 

1. Low-quality qubit (Ti <^ T): phase shift probes the 
parametric inductance of qubit. 

First case which allows to detach the equations for the qubit 
and resonator, is when all the qubit's characteristic times, and 
in particular the relaxation time Ti, are smaller than the tank's 
period T — 2tt/ujo- Then the equations for the tank voltage 
can be averaged over the period of fast oscillations. Then the 
time derivative of the flux (l>e, induced by the qubit in the tank 
circuit can be described as 



(27) 



where $ = $dc + MIi^ is the flux in the qubit's loop, which 
consists of the time-independent part $dc and of the flux, in- 
duced by the current /l in the tank's inductor This can be 
rewritten by introducing the effective inductance of the qubit, 
= dlqh{^)/d^, and the characteristic inductance value 
L = A/2/:-i. Then $e = Z(/l)/l and for the tank voltage 
we have V = Loh - = (-^0 - £(/l))/l- In the first 
approximation in fc^ in the expression $e = L{I\^)I\^ we can 
insert /l found from this equation 



/L(i) j VdtK am{ojpt + S). (28) 

Lo J oJoLo 

Then from Eqs. (I23H24I I we have i?q ^ = (hence R^s = Ro) 
and 



Lo L 



£ ^{vjip) cos'^ ipdip, 



(29) 



Aw = Woff(w) - Wo 



uoLo 
2Lq{v) ■ 



(25) 



The phase shift S and the amplitude v depend on the prob- 
ing frequency detuning ^0 = and the qubit state (via 
Lq and Rq). In the stationary regime they are given by the 
solution of the system of equations 



tan (5 = 2Qo — 



V — IaRcS cos 5, 



(26) 



where the qubit's effective inductance is defined by the total 
flux $, piercing the qubit's loop 



5/qb($) 



(30) 



Then for the phase shift 5 and the voltage amplitude v we 
obtain Q 

tan(5 2(3o^o + QoT^, w ~ /a-Rq cos (5, (31) 
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which is the generalization of the resuh of Ref. f9S'\ for the 
case when the qubit can be in the superpositional state, which 
is taken into account here by the expectation value of the cur- 
rent /qb- If the bias current amplitude 1^ is small enough to 
be ignored in Eq. ( l30l l. then 



c- 



qb 



t&nS « 2Qo£.o + k^QoLC-^, v f^i IaRqCOsS. (32) 

At the resonant frequency — 0, the phase shift 6 is pro- 
portional to the inverse inductance of the qubit C^^. Here it 
is worthwhile to emphasize the expression for the parametric 
inductance, which is expressed via the derivative of the expec- 
tation value of the current in the qubit's loop /qb = — /p (c^). 



° $0 9/dc ■ 



(33) 



2. Higher-quality qubit (Ti < T): parametric resistance 
due to qubit's lagging. 

Another illustrative situation, where the qubit's dynamics 
can be considered separately from the resonator's one, is the 
case when the qubit relaxation time Ti is of the same order as 
the tank's period T, namely, Ti < T. The qubit's response to 
the resonator probing signal can be phenomenologically de- 
scribed by introducing the lagging time t' — t ~ Ti, so that 
instead of Eq. ( |27] | we have 



(34) 



In this way, the qubit's response depends on the current in the 
tank /l = /l (f), which is given by 



{C sm{ujpt + S) - S cosiupt + 6)) (35) 



with S = sin(a;pri) and C ~ cos(wpTi). For the small 
bias current Eqs. (I21H22) and (I34H35) result in the following 
expressions for the parametric inductance and resistance 



i?o/ Rq 



C-k^LC-\ 
-S-k'^QoLC- 



(36) 



By analogy with Eq. ( |33] |. the latter phenomenological equa- 
tion can be rewritten in the form explicitly demonstrating its 
quantum character. 



-rtq — O 



Lip d (0-2 



Ro $0 df, 



(37) 



dc 



By making use of Eq. ([S]), we obtain that the energy trans- 
ferred from qubit into the resonator (or, out of the resonator, 
for the opposite sign) during one period is 



W : 



-nojpv'^Rq ^ . 



(38) 



We emphasize here that both the parametric inductance and 
resistance in Eq. ( |36] | are proportional to the qubit's induc- 
tance C. Then, one obtains equations for 5 and v, which 



are simplified in the first approximation in k'^QQLC~^ . In 
this case for the probing frequency equal to the resonant one, 
^0 — 0, the resulting formulas are 



tanS K C ■ k^QoLC-\ 
-4r « 1 + ^ • k'^QoLC- 



(39) 



Note that both the phase shift and amplitude are related to the 
qubit's effective inductance £, which explains their similar 
behavior in experiment. These equations are useful for the 
analysis of the experimental results, as it will be demonstrated 
in Section 4. 



Capacitive coupling with nanoineclianical resonator. 
Parametric capacitance 

Consider now the charge qubit capacitively coupled to a 
resonator In this case, like in the one considered above, the 
resonator can be the tank circuit. Alternatively, the resonator 
can be a nanomechanical resonator (NR), as in Ref. [i43i1 . For 
the illustrative purpose, we consider here this latter case. 

The split-junction charge qubit (shown in red in Fig. |3]l 
consists of a small island between two Josephson junctions 
(also called Cooper-pair box), whose state is controlled by the 
magnetic flux $ and the gate voltage Vcpb + Vmw- Here 
VcPB is the dc voltage used to tune the energy levels of the 
qubit and Vmw = VJj sin uit is the microwave signal used to 
change the energy-level occupations. The driven Cooper-pair 
box is described in the two-level approximation by the Hamil- 
tonian in the "charge" representation, Eqs. (I13m4t . where 
the tunnel splitting A is equal to the Josephson energy con- 
troUed by the magnetic flux $: A = E'jo |cos(7r$/$o)|- The 
charging energy and the driving amplitude are the follow- 
ing £o = ~8Ec{ng — 1/2) and A — —'HEcn^, where the 
Coulomb energy Eq — e^/2Cs is defined by the total capac- 
itance Cs = 2Cj + CcpB + Cnr and the effective Josephson 
capacitance is introduced 2Cj = Cji + Cj2, the dimension- 
less driving amplitude ~ CcpBV'^/2e. The dimension- 
less polarization charge Tig = rtNR + Jt-cpb is the fractional 
part of the respective polarization charges in two capacitances: 
"-NR = {C'NRVNR/2e} andncpB = {C'cpBV'cpB/2e}. 

The Cooper-pair box here is formed by four capacitors, 
Cji, Cj2, CcpB, and Cnr. One of the plates of the latter 
capacitor is formed by the NR. The displacement of the NR x 
is much smaller than the distance d between the plates. Then 
the capacitance between the NR and the qubit reads 



Cnr^x) « Cnr + 



9C 



NR 



dx 



= C 



NR 



(40) 



Here Cnr stands for the capacitance value at the zero dis- 
placement. The displacement of the NR influences the qubit 
through the changes in the polarization charge; to make it sig- 
nificant, a large dc voltage Vnr is applied. On the other side. 
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GNR RF 



GNR 



the measurable value, the tank circuit phase shift at resonance 
frequency, = 0, 18811 



tan (5 w Qi 



(44) 



cf. Eq. ( I3TI 1. where the phase shift probes the parametric in- 
ductance. In section 4 it will be demonstrated how these ex- 
pressions can be used for the description of the realistic sys- 
tem. 



FIG. 3: (Color online) Charge qubit probed by a nanomechanical 

resonator. The charge qubit is the Cooper-pair box, controlled by the 
magnetic flux $ and the gate voltage Vcpb -|- Vmw- The resonator 
probing the qubit's state here is the NR, which is characterized by the 
displacement at the midpoint x. The voltage-biased NR is measured 
through its resonance frequency shift Aljnr. ISSll 



the NR is biased by dc and rf voltages Vgnr and Vrf through 
the capacitance Cqnr- 

One of the approaches to describe the system qubit- 
resonatoris to introduce the parametric capacitance as follow- 
ing (for more details see Ref. jssll ). Let us introduce the effec- 
tive capacitance, as it is demonstrated in Fig.[Tfc), by differen- 
tiating the charge Qnr of the capacitor Cnr I.99h101i] : Coff = 
dQ^K/dV^B.. Then, for the charge Qnr = {V^k - V[)Cnr 
with the island's voltage given by Vi = 2e{ng — {n))/C^, we 



obtain CcS = Cg, 
capacitance 



Cq, which consists of the parametric 



and the geometric capacitance Cj 



Cnr(Cs — C] 



NR) 



gcom 



2C.tC, 



JL^NR 



2Cj + Cnr 



NR., 



(41) 



(42) 



where the approximations are valid for Ccpb ^ Cj , Cnr 
and Cnr <C Cj respectively. Then one can consider the force 
Fnr , which acts on the NR from the left electrode, as the elec- 
trostatic force from the effective capacitance [see Fig. [TJc)]: 
Fnr = i^csV^^). Then the term with the parametric 
capacitance, in which C^j^ « C^p (1 + x/^)^, results in the 
following resonance frequency shift of the NR 



wnr 



C2 



dn„ 



(3d {a, 
2 dn^ 



(43) 



TOW^pCs 



CnrVnr, \ ' 



We would like to note that the results obtained for the sys- 
tem qubit-NR can be definitely extended to other systems. For 
example the charge qubit can be coupled to a tank circuit in- 
stead of a NR. In contrast to the inductive coupling, consid- 
ered in the previous subsection, here we mean capacitive cou- 
pling. Then it is straightforward to obtain the expression for 



3. DYNAMICAL BEHAVIOR OF A TWO-LEVEL SYSTEM 

Application of the semiclassical theory, presented in the 
previous subsection, to the description of the qubits-resonator 
system makes possible to separate the slow dynamics of the 
resonator from the fast dynamics of the qubits system. This 
allows to consider first the dynamics of a qubit or a system of 
qubits. Then, the resonator can monitor the state of the system 
of qubits. In this section we will outline the description of the 
multiphoton processes in a qubit, while the presentation of the 
specific results is the subject of the next two sections. 

Initialization and manipulation of the qubit's systems re- 
quire certain external signals. The principal features of the 
driven system are captured for the harmonic driving, Eq. ([T4l i. 
to which we limit our consideration. Different theoretical ap- 
proaches can be used for a driven tw o-level sys tem, which is 
described in the books and reviews 1*102- 107 1. The choice 
of the formalism depends on the formulation of a problem 
and on the parameters of the system, such as the bias offset 
Eq, driving amplitude A and frequency lu. The clear descrip- 
tion can be given for the temporal dynamics in the so-called 
adiabatic-impulse model, where the driven evolution is con- 
sidered adiabatic far from the avoided-level crossings with the 
impulse-type Landau -Zener tran sitions, when the energy dis- 
tance is minimal ll37[ll08Lll09ll . As the result of this theory, 
the overall dynamics is described by the long-time Rabi-type 
oscillations of the level occupation probabilities with the step- 
like features due to the Landau-Zener transitions. 

Another technique, which can be more convenient for the 
resonant d riving, is the rotating wave approximation (RWA) 
1 1 lOl - 1 12 1. It consists in neglecting the rapidly oscillating 
(non-resonant) terms. The common approach for making use 
of this approximation is taking small driving amplitudes, A <C 
AE. Then, the first-order consideration gives usual Rabi os- 
cillations of the level occupation probabilities close to the po- 
sition of the one-photon resonance, where w « AE/h. In 
the fc-th approximation, the resonant excitation appears close 
to the parameters, where the energy of k photons matches the 
qubit's energy distance HI 01 



khuj = AE. 



(45) 



The time evolution is described by the multiphoton Rabi oscil- 
lations 1 113], while the time-averaged upper-level occupation 
probability has the Lorentzian shape with the maximum equal 
to 1/2 at the exact resonance defined by Eq. ( |45] |. 
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W ith increasing the driving amplitude the resonances shift 
Ol 1411 from their positions given by the perturbation theory and 
defined by the exact multiphoton relation ( |45] l. The first-order 
correction to the position of the resonances is the so-called 
Bloch-Siegert shift [IJ; itwas demonstrated for the supercon- 
ducting qubits in Ref. Ill 1511 . Thus, in general, the position of 
the multiphoton resonances is amplitude-dependent. 

For the description of the strongly driven qubits, another 
formulation of the RWA can be used. There, the mini- 
mal energy level splitting A is th e small par ameter, namely, 
it is assumed A < \fAbj Q [TTel [iItIi . Then the fc- 
photon excitation appears close to the resonant parameters, 
given by the relation £o = Khuj. There, the upper-level 
occupation probability Pup(^) oscillates with the frequency 
51r — ^ (eo — khoj)^ + Al with the renormalized splitting 
Afc = AJk (A/uj); Jk is the Bessel function. The time- 
averaged probability in the vicinity of the fc-th resonance is 
given by 



1 



2 (eo - khuoY + Al ■ 



(46) 



Being time averaged, the Rabi oscillations are described by 
the Lorentzian dependence of the upper-level occupation on 
the s ystem's parameters (the bias or the driving frequency) 
01 1811 . Here arises an interesting and important problem of 
distinction of the respective quantum oscillations from their 
classical counterparts, which are the^parametric resonances. 
This was the subject of Refs. ilT9Ml2lll . 

The most straightforward approach for the numerical de- 
scription of the dynamics of a two-level system is the solu- 
tion of the Schrodinger equation rfl22ll . Then, the influence of 
the dissipation can be taken into account phenomenologically 
by introducing energy and phase relaxation times, Ti and T2, 
and solving the respective Bloch equation lll02ll . Instead, in 
the more general approach, the dissipative environment can 
be described as an ensemble of oscillators, which would re- 
sult in the Bloch -Redfield equation for the reduced density 
matrix 1 123 , 124^. This latter formalism will be demonstrated 
in Sec. 5 being applied to the specific case of the two-qubit 
system. 

Note that the multiphoton transitions can also be driven 
by the bichromatic field, when the energy level distance AE 
is matched by the energy of several photons of one (say, 
microwave-) frequency plus several photons of another (say, 
radio-) frequency. Such transitions were studied both in mi- 
crosco pic system s [T25il . and in the Josephson-j unction 
qubits 1261 - 1 12811 . Also for the case of a flux qubit it was 
demonstrated that the persistence of Rabi oscillations can be 
supporte d by e ither the low-frequency signal lll29ll or induced 
by noise lfl30ll . 



4. EXCITATION OF A SUPERCONDUCTING QUBIT 

Let us get back to the qubit-resonator systems. In the pre- 
vious Section we have discussed a modification of the qubit 



states (and therefore its observables) under different types of 
excitations. A natural next step is to analyze the correspond- 
ing (via qubits) change of the resonator properties. In this 
section we demonstrate this by presenting respective theoreti- 
cal results for different realizations of the qubit-resonator sys- 
tems, making use of the theory presented in the previous two 
sections. The emphasis is made on demonstrating the consis- 
tency of the theoretical results with the experimental ones. 



Inductance of superconducting qubits 

Consider a qubit biased with a DC flux $dc and driven with 
an AC flux <I>ac sinwt, introducing /dc — ^dc/^o ~ 1/2 and 
/ac = <&ac/<&o- In order to get the effective inductance C, 
as defined by Eq. (l30l l. we have to calculate the average cur- 
rent in the qubit: /qb — (/) = Tr (pi), where / = /pCT^ 
is the current operator defined with the amplitude Ip and the 
Pauli matrix a^- We calculate th e red uced density matrix 
p with the Bloch equations [102l 112211 which include phe- 
nomenological relaxation times, Ti and T2. It is convenient 
to express the density matrix in the energy representation: 
p — (1/2) (to + + YTy + Ztz), where are the Pauli 
matrices for this basis and tq stands for the unity matrix. The 
value Z ~ {%) is equal to the difference between the popula- 
tions of the ground and excited states. 

Let us find now the explicit expressions for the effec- 
tive qubit's inductance f or both the interferometer-type (split- 
junction) charge qubit 1113 ll 113211 and flux qubit [94]. For 
the interferometer-type charge qubit, as considered in detail 
in Ref. fsd], the circulating current Iq is flux-dependent and 
Eqs. ( |32] | show that there are two terms contributing in the 
tank circuit's phase shift, 

fc2QL f dlo 



tan 6 



$0 \dU 



-Z + Io 



dZ 
dfdc 



(47) 



In a classical system (where the current has a definite direc- 
tion) or in the ground state, the difference between the energy 
level's populations is constant, Z = const, and the second 
term in Eq. ( |47] i is zero. In contrast, for the quantum sys- 
tem the interplay between these two terms is essential. At 
this point it is worthwhile to notice that the second term can 
dominate at resonant excitation, as it was the case in the work 
[30] (see also below). This means that the second ("quan- 
tum") term can significantly increase the sensitivity of the 
impedance measurement technique, as compared to the clas- 
sical situation described by the first term in Eq. (l4Tt . 

Consider now the case of a flux qubit. The current operator 



is defined in the flux basis [94l|, / = /pcr^, where Ip stands 
for the amplitude value of the persistent current, and hence 
the value (cr^) defines the difference between the probabilities 
of the clockwise and counter-clockwise current directions in 
the loop: (ct^) = - = 2P; - 1. Then with Eqs. ^ we 
obtain 



(48) 



9 



In the energy representation we rewrite Eq. (148) 

tan<5«fcQ— ^(^— (49) 



Here AE = + (/p$o/dc)^ is the distance between the 

stationary energy levels. 

After the time-averaging over the driving period 2tt/lu, this 
expression is written as following 



tan S : 



Lll ( A3 



A \l\E'^ 



(50) 



If a qubit is resonantly excited with the driving frequency ui, 
then the partial energy levels occupation probability Z has 
the Lorentzian-shape dependence on /dc- It follows that the 
derivative dZ/ dfdc takes the shape of a hyperbolic-like struc- 
ture, i.e. it changes from a peak to a dip in the point of the 
resonance at AE{fdc) ~ fcfiw. 



Equilibrium-state measurement 

For the description of the measurement of a flux qubit in 
the thermal equilibrium one has to put X = and Z = 
tanh {AE/2kBT) in Eq. (US, 



tan (5 ; 



A3 



2 



(- 



A d 
AE^^'dUc 



tanh 



AE 



2kBT 
(51) 

The ground-state measurement at fee T <^ AE is described 
with X = and = 1, which means replacing the hyper- 
bolic tangent in Eq. JSTl ) with the unity. The formula ( ISTT i 
for the ground state obtained by differentiating the probabil- 
ity Pi, Eqs. ( l48H49t . coinci des w ith the earlier obtained re- 
sults (see Eqs. (3-4) in Ref. il33ll ). The resulting tank phase 
shift is shown in Fig. |4] for the following parameters taken 
from Ref. ijBl: A//i = 1.3 GHz, Ip<Pa/h = 930 GHz, 
uJo/27T = 32.675 MHz, L/p/$o = 0.0055, M/L = 0.725, 
Qo = 725, k = 0.02. 

The accurate account of Z in Eq. ( BTT i allows to describe 
both the suppression and widening of the zero-bias dip (that 
is at frir = 0) as it was experimentally demonstrated in 



Ref. 1113411 . Indeed, the suppression of the zero-bias dip (at 
/dc = 0) is described by the first term in Eq. (ISTi . The widen- 
ing is due to the second term that comes from differentiating 
the hyperbolic tangent; this term becomes relevant for temper- 
atures larger than A, and results in the exponential rise of the 
width for T > T* — A/k-Q, as demonstrated in the inset in 

Fig. a 



Resonant transitions in the charge qubit 



In Ref. [301] the resonant excitation of the interferometer- 
type (split-junction) charge qubit was demonstrated experi- 
mentally and described theoretically. In accordance with the 
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FIG. 4: (Color online) The equilibrium-state measurement. The 

dependence of the tank phase shift on the flux detuning /dc = 
^dc/'^Q — 1/2, when the qubit is thermally excited. The curves are 
plotted for kBT/h = 0.2, 0.5, 0.7, 1, 2, 4, and 8 GHz. Left inset: 
corresponding experimental results ||134|| . Right inset: temperature 
dependence of the width A/dc of the dip at half-depth in the phase 
shift, shown in the main panel. Ii87i1 



formula (|47] | one expects the resonances to appear differently 
when either first or the second term is dominated. To demon- 
strate this, in Fig.|5]we plot the dependence of the tank circuit 
phase shift 5 both as the function of the dimensionless bias 
voltage Tig = CgV^/2e and of the dimensionless magnetic 
flux detuning /dc. For the former case the value /dc — was 
taken, where /q 0. This results in disappearance of the sec- 
ond term in Eq. (l47T i. and the resonant excitation of the qubit 
is visualized with the Lorentzian peaks in Fig. |5la,b). When 
the second term is dominant, the multiphoton transitions in the 
qubit result in the peak-and-dip structures in the dependence 
of the phase shift 5 on the flux, Fig.|5jc,d). 

Theoretical fitting of the experimental graphs, as for exam- 
ple shown in Fig.|5] allows for defining the qubit's parameters, 
which is the multiphoton spectroscopy. The parameters found 
were the following: the Josephson energies for the two junc- 
tions Eii/h ~ 40 GHz and E^/h ~ 34.5 GHz, the island's 
Coulomb energy Eq //i ~ 5 GHz; the relaxation and decoher- 
ence rates T,^x^^/{Ec/h) = 0.03 and T^/{Ec/h) = 0.05, 
which correspond to the following relaxation and decoherence 



times: T, 



rcla 



r 



rcla 



7 ns and Tj, = T', 



4 ns. 



Figure |5] also demonstrates how the position of the reso- 
nances depend on the driving frequency lj and how the mul- 
tiphoton resonances appear with increasing the driving power 
Uac- Namely, first, in Fig.|5ja,b) the varied parameter is the 
frequency a;/27r, which from the bottom to top curves is 6.5, 
7.1, 8.1, and 9.1 GHz; the driving power is the same for all fig- 
ures nac — 0.3 and the flux was fixed at S = n. And, second, 
in Fig. |3c,d) the curves correspond to the varied parameter 
driving power: in experiment being power of excitation (from 
bottom to top: —80, —60, —57 dB) and in theory being am- 
plitude riac (from bottom to top: 0.1, 0.2, 0.4); the frequency 
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FIG. 5: (Color online) Resonant excitation of the charge qubit 
probed by the tank circuit. The phase shift 5 of the tank circuit 
coupled to the charge qubit, calculated theoretically (left) and mea- 
sured (right). Panels (a) and (b) show the dependence on the gate 
voltage, while in (c) and (d) the dependence on the flux is demon- 
strated. Black and gray arrows in (c) demonstrate the positions of 1- 
and 2-photon resonant transitions, and the arrows in (d) mark 1-, 2-, 
and 3-photon excitations, fsoll 



there was fixed, a;/27r — 7 GHz. 



One- and multiphoton transitions in the flux qubit 

As we have seen in Section 3, both the tank voltage phase 
shift S and amplitude v can be used to monitor the resonant ex- 
citation of a superconducting qubit. In Fig.|5]we demonstrated 
this with the observation of the phase shift S of the tank cir- 
cuit coupled to the charge qubit. Now we consider one- and 
multi -photon resonant excitations of a flux qubit, and the non- 
monotonic dependence of the tank voltage amplitude v will 
visualize the resonant transitions in the qubit. 

Consider first the spectroscopical measurement, where 
the flux qubit is driven with the low-amplitude AC flux. 
We expect resonant excitation of the qubit when the driv- 
ing frequency matches the qubit's energy difference, hw — 
AE{fdc)- In the experimental case the positions of these res- 
onances at a given driving frequency allow to determine the 
energy structure of the measured qubit [41 1. 

In Fig. |3b,c) we demonstrate the dependence of the tank 
voltage amplitude v on the bias flux /dc at Up — uto for differ- 
ent driving frequencies: a;/27r = 3.5, 5, and 18 GHz, which 
is explained by the energy diagram in|6ja). The results of the 
related experiment, Ref. |41], are presented in Fig. |6lc). The 
parameters for calculations were taken as following: the tun- 
neling amplitude A/h — 3.5 GHz, the energy bias Ip^o/h = 
700 GHz, the temperature /cbT/Zi =1.4 GHz, the relaxation 
rate Ti/h = 0.7 GHz, the dephasing rate T2/h = 0.7 GHz, 
and the value which describes the coupling between the qubit 
and the tank circuit k'^Qo{LIp/^Q) — 2.6 • lO^'^. The curves 
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FIG. 6: (Color online) Low-amplitude one-photon resonant exci- 
tation of a flux qubit. (a) Energy levels E±{fdc) matched by the 
driving at frequencies shown by the numbers and the arrows of the 
respective length, (b) and (c) Theoretically calculated and experi- 
mentally measured amplitude of the tank voltage v versus flux de- 
tuning /dc for different driving frequencies. (The upper curves are 
shifted vertically.) The one-photon excitations at ci;/27r = 18, 5, and 
3.5 GHz, demonstrated in (b) and (c), are explained by the arrows 
to the left in the energy diagram (a), while the arrows to the right of 
the len gth aj/2 7r = 4.15 GHz explain the multiphoton resonances in 

Fig.iniiiii 



were plotted for the driving amplitudes /ac • 10"^ = 1, 1.5, 
and 3 from bottom to top. The phenomenological lagging pa- 
rameter was taken S — 0.8. Figure |6] demonstrates the effect 
described in section 3: for S both the phase shift S and 
the amplitude v depend on the qubit's inductance C^^, which 
results in the alternation of peak and dip around the location 
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Note that for the lagging parameter close to 1 (here S = 
0.8) the changes in the phase shift in Fig. [Tja) are small at 
the resonance frequency (along the dotted line at Up = ujq) 
while the voltage amplitude in Fig. |7jb) changes substan- 
tially, see formulas ( |39] l. And this is actually demonstrated 
in Fig. |6lb,c). Such changes of the tank effective resistance 
or, equivalently, quality factor were studied in Ref. [|90ll for 
the fully quantum-mechanical model of the qubit-resonator 
system. We note that this can be alternatively described with 
the semiclassical model, presented here. This model gives re- 
sults consistent with the experimental ones, e.g. Figs. |6] and 
|7] which imply the energy transfer between the qubit and res- 
onator according to Eq. (l38T l. More details about this energy 



-0.02 -om f 0.01 0.02 -0.02 -0.01 , 0.01 0.02 



transfer, known as the Sis yphus damping and amplification, 
can be found in Refs. Q, (u^ . 



FIG. 7: (Color online) Multiphoton excitations of a flux qubit. 

Theoretically calculated dependence of the phase shift S (a) and the 
amplitude v (b) on the bias current frequency ojp and the flux detun- 
ing /dc- (c,d) Experimentally measured phase shift 5 and the ampli- 
tude V. L87J 



of the resonances. 

In Fig. |3a,b) we present the calculated phase shift S and 
the amplitude v as functions of the probe current frequency 
ojp and the flux detuning /dc with the phenomenological lag- 
ging parameter S for the strongly-driven flux qubit with the 
parameters being the same as for Fig. |6]and with the values 
for the driving amplitude and frequency: /ac = 8 • 10^"^ and 
w/27r = 4.15 GHz. The top panel presents theoretical cal- 
culations, which is in good agreement with the experimen- 
tal observations, presented in the bottom panel. Fig. |2lc,d). 
The dashed white line shows the tank resonance frequency 
ujp/2TT = u!o/2tt = 20.8 MHz. The positions of the mul- 
tiphoton resonances is explained by the arrows to the right in 
the energy diagram, Fig.|6la), at AE{fdc) — khw with k — 1, 
2, 3, and 4. 



Then, in Fig. |8]we present the dependence of the tank volt- 
age phase shift 5 on the microwave amplitude /ac and the DC 
flux bias /dc. This double quasi-periodical dependence (on 
both the energy bias and the driving amplitude) is called the 
Landau-Zener-Stuckelberg interferogram iIstIi . The parame- 
ters were taken the same as for Fig. i] and uj/2tt = 4.15 
GHz. The left panel in Fig. [8] presents the theoretical in- 
terferogram from Ref. 11871] while the right panel is the ex- 
perimental one, Ref. 14111 . In Fig. |8]the multiphoton reso- 
nances at discrete DC bias /dc (which controls the distance 
between energy levels) are clearly visible. These resonances 
appear when the energy of k photons matches the qubit's en- 
ergy levels, khw « AE{fdc)- The quasi-periodical character 
of the dependence on the AC flux amplitude /ac is known as 
Stiickelberg oscillations. The comparison of such graph to the 
experimental analogue allows the relation of the microwave 
power to the AC flux amplitude /ac to be determined, which 
is the calibration of the power For this, either the estimation 
of the period of Stiickelberg oscillations, shown by the black 
arrow, or adjusting the interference pattern slope, shown by 
the white hne, can be used. 




FIG. 8: (Color online) LZS interferometry for the flux qubit probed by the tank circuit. The calculated (a) and measured (b) dependence 
of the tank phase shift on the flux detuning /dc and on the driving flux amplitude /ac. L4L.87i1 
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Interferometry with nanoresonator 

The formalism developed in Sec. 3.3 allows to describe 
the system of the driven qubit coupled to the NR. As it was 
demonstrated in Ref. ifssll . two different approaches, called 
direct and inverse LZS interferometry, are of interest. In the 
direct interferometry the qubit state is probed via the NR's fre- 
quency shift, as in Ref. i43ll . while in the inverse interferome- 
try the impact of the NR's state on the qubit's Hamiltonian is 
studied. 



The direct LZS interferometry was calculated in Ref. 
as the resonator's frequency shift Awnr versus the energy 
bias rig and the driving amplitude riu- The agreement with 



the experimental result of Ref. 114311 demonstrated that the 
semiclassical formalism is valid for a description of the mea- 
surable quantities. In Ref. 188,1 it is also demonstrated how 
the analogous interferogram can be calculated for the qubit- 
tank circuit system in relation to the experiment of Ref. ll39ll . 
Such a description allows to correctly find the position of the 
resonance peaks in the interferogram and to demonstrate the 
sign-changing behaviour of the parametric capacitance, which 
relates to the measurable quantities. 

For the formulation of the inverse problem, let us consider 
the qubit's bias Eq as a function of the NR's displacement x. 
For small x ^ ^ we have the expansion ( |40] |. which results in 
the decomposition of the bias eo(a;) sa e5("-g)+^£o(a;), where 
e*(ng) = 8Ec (rig - 1/2) and Seo{x) = SEcmnx/C The 
Hamiltonian of the qubit ( fT3] l with the parameter-dependent 
bias Eg {x) allows to consider the following problem. Let us 
assume that the qubit's state (its wave function, upper level 
occupation probability, Rabi frequency, etc.) is known (i.e. 
this is measured by a device, which we do not consider here 
for simplicity). Given the known qubit's state, the aim is to 
find the Hamiltonian's parameters. Particularly interesting is 
the parameter-dependent bias eo{x), which can give the infor- 
mation about the position and amplitude of the oscillations of 
the NR. 

And now, in the general co ntext, the "reverse engineering" 
problem in the spirit of Refs. lll36[ll37ll can be studied, where 
one is interested in finding the driving Hamiltonian for a given 
(desired) final state. On the other hand, in Ref. |88] the au- 
thors provide the basis for measuring the NR's position x by 
means of probing the qubit's state, while x — x{t) is con- 
sidered a slow time-dependent function. There, the empha- 
sis was made on finding optimal driving and controlled offset 
(Eq) parameters for the resolution of the small bias component 
6eq. It was assumed that the dynamics of the parameter x is 
slow enough not to be considered during either certain period 
of the qubit's evolution or even during the setting the station- 
ary qubit's state. The aim was to find a sensitive probe for 
small 6eq. As the ultimate sensitivity, the essential changes 
of the qubit's state for small changes of Seq were required. 
The problem, formulated in this way, was solved in Ref. ll88ll 
for different illustrative driving regimes: one-, double-, and 
multiple-passage regimes. 




Bath of harmonic oscillators 



FIG. 9: (Color online) Scheme of two coupled qubits. The two flux 
qubits are coupled to each other, to the dc and /iw lines, as well as to 
an unavoidable dissipative environment. The convenient model for 
description of the environment is the bath of harmonic oscillators. 
The system of two coupled qubits is also assumed to be coupled to 
the measuring resonant circuit (which is not shown here), as in Fig.|2] 



5. MULTI-QUBIT SYSTEMS 

Equations for a system of coupled qubits 

The effective Hamiltonian of the system of n coupled flux 
qubits is 



H 



i=l 



A,; 



J, 



i,3 



(52) 



where Jij is the coupling energy between qubits, and <Tx\ 
a't^ are the Pauli matrices in the basis {\\.) , W)} of the cur- 
rent operator in the z-th qubit. The current operator is given 
by: li = — Jp tJz , with Ip the absolute value of the per- 
sistent current in the i-th qubit; then the eigenstates of cr^ 
correspond to the clockwise (ctz = — and counter- 
clockwise (az It) = It)) current in the i-th qubit. The tun- 
neling amplitudes are assumed to be constants. The biases 
Ei = 2Ip^ ^of^''^ {t) are controlled by the dimensionless mag- 
netic fluxes /'^'■\t) = <I>i/<i>o - 1/2 through i-th qubit. These 
fluxes consist of three components, 

f^'Ht) = f, + ^ + Usmujt. (53) 

Here fi is the adiabatically changing magnetic flux, experi- 
mentally applied by the coil and additional DC lines. The sec- 
ond term describes the flux induced by the current /l in the 
tank coil, to which the i-th qubit is coupled with the mutual in- 
ductance Mi. And /ac sinwi is the harmonic time-dependent 
component driving the qubit, typically applied by an on-chip 
microwave antenna. Equation (|52] | can be reduced to the two- 
qubit system. This system is shown in Fig.|9] 

To describe the two-qubit system, it is convenient to present 
the density matrix in the following form 



RaO 



P = 



I (To 



Rap 
Rob 



^00 ^ 
1 cr/3 = — ;— CTq (8" CTo 



-CTo 8) (7b 



4 
Rab 



-<Jn 



(54) 
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which was shown to be suitable for both the definition and 
the calculation of the entan glement a nd other characteristics 
in multi-qubit system, e.g. 1 138 , 139ll . Here a, /3 = 0, x, y, z 
and a, 6 = x^y^z; the summation over twice repeating in- 
dices is assumed. The two vectors i?ao and R^h, so-called 
coherence vectors or Bloch vectors, determine the properties 
of the individual qubits, while the tensor Rab (the correlation 
tensor) accounts for the correlations. 

The important characteristic of the state of the coupled- 
qubits system is its entanglement. There are different ap- 
proaches to the quantification of the entanglement 1 140]. One 
of th e often used possibilities is the so-called concurrence 
1 14111 . Another convenient for calculations approa ch is to in- 
troduce the measure of entanglement as following 11138 1 



£ - ixr (M^M) , Afab = Rah - RaoRob- (55) 

This entanglement measure fulfills certain requirements, in 
particular, £ = Q for any product state and £ = \ for any 
pure state with vanishing Bloch vectors Rao and i?oh, corre- 
sponding to maximum entangled states. 

To describe dynamics of the density matrix we will first dis- 
regard the relaxation processes. This can be described by the 
Liouville equation, ihp = [H, p], which is generally speaking 
a complex equation. To deal with the Liouville equation, it is 
convenient to use the parametrization of the density matrix as 
described by Eq. (|54] |. Due to the hermiticity and normaliza- 
tion of the density matrix, Rap are real numbers and Rqq = 1. 
Then the Liouville equation can be written in the form of the 
system of 15 equations for i?Q,^ |541] 
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+ SjS^SniJRnO + Si3^3njJRon, 



(56) 



where B^'^ — (— A^, 0, — e^) and Smni is the Levi-Civita sym- 
bol. 

Consider now the measurable value, which is the res- 
onator's phase shift. As we discussed in Sec. 2, it relates to the 
effective inductance of qubits system. The formula obtained 
for singl e qubits can be generalized for the two-qubit system 
119711 14211 . Then for the case of low-quality qubits, when their 
characteristic times are smaller than the tank's period, at the 
resonance frequency (^o = 0), expression for the phase shift S 
in terms of the parametric inductances Lq ■* can be written as 
following 



tan 5 



•3° 7W' 

1=1,2 



(57) 



^0 _ ,2^i /'-I _ ( ^ 



d 



\d^a d^bj '''' 



In what follows this expression will be used to calculate the 
phase shift S, which maps the qubits' state. 




/a 



8 [rad] 



FIG. 10: (Color online) Spectroscopy of the two-qubit system. The 

measured dependence of the phase shift S on the flux biases fa and 
/{,: (a) ground-state measurement (without microwave excitation). 
(b,c,d) with weak microwave excitation at the driving frequencies 
uj/2tt = 14.1, 17.6 and 20.7 GHz. (H 



Weak-driving spectroscopy 

In Sec. 4 we have considered how the measurements of the 
single qubits allow to determine their parameters: the tun- 
neling amplitudes A and the persistent currents /p. It was 



demonstrated 114 111 that for defining the parameters of single 
and multiple-qubits systems both the ground-state measure- 
ments and excited-state spectroscopy can be used; the consis- 
tency of the results of the two approaches was shown. Now we 
will demonstrate this for the case of the system of two coupled 
flux qubits described by the Hamiltonian (|52] |. 

First, the one-qubit parameters are defined. For this, sup- 
pose qubit a is the one biased far from its degeneracy point in 
such a way that is large in comparison with the other en- 
ergy variables. Then, qubit a has a well defined ground state 



with averaged spin variables {^^z^^^ = 1 ™d i^'^x"' j = 
which can be averaged out of the two-qubit Hamiltonian i 

reducing it to: iJ2qbs,rcd = -Abcri^^ /2 ~ {sb 
Apart from the offset in the bias term due to the coupling, this 
is identical to the single qubit Hamiltonian. This offset can 
be easily compensated and measured allowing the determi- 
nation of the coupling energy J lll42ll . The qubit parameters, 
A;, and Ip"^ , are determined from either the ground-state mea- 
surement or the excited-state spectroscopy, as it is described in 
Sec. 4. Analogously, biasing qubit b far from the degeneracy 
point the parameters for qubit a, Aa and Ip"'\ can be deter- 
mined. Next, the coupling energy J was determined from the 
offset of the qubit dips from the ^a/b = lines, visible in the 
pure ground-state measurements presented in Fig. [TOla). 

Then the qubits were driven by magnetic fluxes <I>ac sin ojt 
with weak driving amplitudes and various driving frequen- 
cies. There, we expect the position of the resonant transitions 
from energy level Ej to an overlying level Ei determined by 
the one-photon relation: AEij{^a, ^b) ~ which appears 
when the distance between the energy levels AEij = Ei — Ej 
is matched by the photon energy hu. In Fig. fTOf b) a frequency 
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in-between both qubit gaps (A;, < Hlu < Aa) was used and 
therefore only the transitions to the first excited state are visi- 
ble. For higher frequencies, also the second and third excited 
states become visible as can be seen in subfigures (c) and (d). 
The theoretically calculated contour lines are superposed in 
Fig. [TOl b-d) for three different frequencies for which the con- 
dition, AEij{^a,^b) = is fulfilled; the energy levels 
Ei — Ei{^a,^b) were found by diagonilizing the Hamil- 
tonian. From the fitting procedure the following parameters 
were found: the tunneling amplitudes Aa{b)/h — 15.8(3.5) 

GHz, the energy biases I^^^^^a/h = 375(700) GHz [7°^^' = 
120(225) nA], the inter-qubit coupling J/h = 3.8 GHz, and 
the value which describes the coupling between the qubits 



and the tank circuit ^ 

A:fgo(L,4'V*o). 
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Direct and ladder-type multiphoton transitions 



FIG. 1 1 : (Color online) Characterizing strongly-driven two-qubit 
system. Calculated and plotted as functions of the bias /dc are four 
energy levels (a), total probability of the currents in two qubits to 
flow clockwise Z (b), the tank circuit voltage phase shift 5 (c), the 
entanglement measure £ (d). (stIi 



We now consider the multiphoton excitations of a system of 
two strongly driven coupled flux qubits. We will describe the 
effects of resonant excitation in the system in terms of its en- 
ergy structure, entanglement measure, and the observable tank 
circuit phase shift. Then we will present results for the multi- 
photon excitation of two types: direct (when multiple photon 
energy khuj matches the energy level difference AEij) and 
ladder-type (when the transition happens via an intermediate 
level). We will demonstrate how this can be used for creating 
the inverse population in the dissipative two-qubit system. 

To describe the system of two qubits subjected to the strong 
driving, the following values were calculated: the energy lev- 
els (by diagonalizing the stationary Hamiltonian), the density 
matrix p (by solving the Liouville equation), the observable 
tank circuit phase shift & (which is defined with the effective 
inductance of the qubits), and the entanglement measure £ by 
making use of Eqs. ( l55H57b . In this way graphs in Fig.fTTIwere 
calculated for th e set of parameters of the two-qubit system 
realized in Ref. |[lil: Aa/h = 1.2 GHz, Ab/h = 0.9 GHz, 
4°'''^$o/^ = 990 GHz, J/h = 0.84 GHz, ^a,b = 1-8 • lO'^, 
and the driving frequency was taken lli/2tt = 4 GHz; also 
the change of the DC flux here was assumed symmetrical: 
fa = fb^ /dc- For simplicity here the relaxation processes 
were ignored (and we will pay special attention to this below) 
and we consider the case when the characteristic measurement 
time Tp = 2tt/ cjp is larger than the characteristic times of the 
dynamics of the qubit. Then the tank circuit actually probes 
the incoherent mixture of qubit's states and the time-averaged 
values of phase shift and entanglement should be considered. 

When the energy of k photons ikhw) matches the energy 
difference between any two levels Ej and Ei, the resonant 
excitation to the upper level is expected. Respectively, the 
arrows of the length 4, 8, and 12 GHz show the places of 
possible one-, two-, and three-photon excitations. The time- 
averaged total probability of the currents in two qubits to flow 
clockwise, Z = i?03 + i?3o, is shown in Fig. niTb) to expe- 



rience resonant excitation. The resonances appear as peak- 
and-dip structures in the phase shift dependence in Fig. fTTT c). 
The time-averaged entanglement measure f in a resonance in- 
creases due to the formation of the superposition of states. 
Fig. [TT? d); this provides a method to control and probe the 
entanglement. 

The experimental study of the strongly driven system of two 
coupled flux qubits is presented in Fig. [12] The left panel is 
the measured voltage amplitude of the tank as a function of 
qubit biases /a and /b. The driving frequencies from top to 
bottom were lo/2tt — 17.6, 7.0, and 4.1 GHz. The multipho- 
ton resonances at AEij {fa, fb) ~ khuj are visualized with the 
ridge-trough lines. We note that the resonance ridge-trough 
lines are disturbed with increasing or decreasing the signal; 
some of these changes are shown with white circles. This 
means changing the effective Josephson inductance in these 
points. The experimental results can be clearly understood 
by comparing them with the energy contour lines, calculated 
by diagonalizing Hamiltonian ( |52] | and presented in the right 
panel of the figure. There, numbers fc — j next to the lines 
mean that the line relates to the energy difference Ej — Ek ■ 

Consider now these multiphoton features in more details. 
In Fig. [TSlb) the black and red lines show the positions of 
the expected resonant excitations from the ground state to the 
first and to the second excited states respectively; the blue and 
orange lines are the contour lines for the possible excitations 
from the first and from the second excited state to the third 
excited state. In Fig.[T3ja) the energy levels are plotted at the 
fixed value of the bias flux through qubit a, /a, as a function 
of the bias flux through qubit b, /b. The arrows are intro- 
duced to match the energy levels with the driving frequencies 
uj/2tt — 17.6 GHz and 7.0 GHz . The black and red arrows in 
both Fig. [T3l a) and Fig.[T2|b) show the position of one-photon 
transitions to the first and the second excited levels. The dou- 
ble green and blue arrows in Fig. [T2]show the position of the 
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FIG. 12: (Color online) Imaging the multiphoton transitions in 
the two-qubit system. The resonant excitation of the qubits system 
is visualized by the tank voltage amplitude, (a, c, e). The position 
of the resonant transitions can be understood by comparing with the 
respective energy contour lines (b, d, f). jssll 



two-photon processes, where the excitation by the first photon 
creates the population of the first and the second levels and 
the second photon excites the system to the upper level. These 
two-photon excitations happen via intermediate levels; com- 
pare the position of these expected resonances in Fig. [TSl b) 
shown with the blue circle and green square. The orange tri- 
angle in Fig.[T2|b) points the ladder- type three-photon excita- 
tion, with one photon to the first excited level and then with 
two photons to the upper level. 

Analogous considerations allow to see in Fig. fTSl c) and (d) 
one- and two-photon resonant excitations to the first excited 
level for the driving frequency lli/2'k = 7 GHz. The two- 
photon resonant excitation is direct and happen without any 
intermediate level. The higher level excitations via the first 
excited state appear due to three- and four-photon excitations, 
as shown with orange triangles and pink asterisk. In Fig. [T2l e) 
the response of the two-qubit system at w/27r = 4.1 GHz 
exhibits 1- to 4-photon excitations to the first excited state, 
which can be recognized by comparing with the black lines 
in Fig. [T2l f). Numerous upper level excitations via the first 
excited level appear as the changes of the signal along these 
lines. 

The transition rates can be quantified by the absolute value 
of the matrix element of the perturbation between the states 
and \E„) 



Turn = \{E„\v\E„ 



(b) V P ^ 



divided by the factor /p^'$o/ac- The transition matrix el- 
ements in Fig. [TSl'b) explain the ladder-type excitations in 
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FIG. 13: (Color online) Ladder-type transitions in the two-qubit 
system. Calculated as functions of the flux /b (at /a = 0.015): the 
energy levels (a), transition matrix elements r„m (b), the occupation 
probabilities Pi (c). |53] 



Fig. fl^ b). Two points, marked by the vertical dashed green 
and blue lines in Fig. [13] describe respectively two interest- 
ing situations. To the right (see along the blue line) the tran- 
sition element between the higher two levels (£'2 and £"3) 
is smaller than between the lower two levels {Eq and E2), 
Tq2 » T23 » T03. In contrast, to the left (see along the green 
line) the transition element between the higher two levels {Ei 
and £3) is larger than between the lower two levels (Eq and 
El), Ti3 > Toi > T03. In both cases the probabihty of 
the direct excitation to the highest level is very small, which 
means that the transitions are induced due to the ladder-type 
mechanism. 

The ladder-type transitions and the population inversion can 
be also illustrated by calculating the energy level occupation 
probabilities by solving the Bloch-Redfield equation (see the 
next subsection for more details); figure [T3l c) was calculated 
with the driving frequency w/27r — 17.6 GHz and ampli- 
tude /ac = 4 X 10^"^. First, the ladder-type resonant exci- 
tation takes place to the left, where the upper level occupa- 
tion probability P3 is of the same order as the intermediate 
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level occupation probability Pi. Second, the inverse popu- 
lation appears to the right, where the upper level occupation 
probability Pi is l arger than the ground state probability Pq, 
see also Refs. li20il42i[l44lll45ll for the study of the popula- 
tion inversion in the systems with single Josephson-j unction 
qubits. These two phenomena are similar to those which ex- 
hibit atoms in the laser field [146]. Furthermore, the expec- 
tation value of the current in i-th qubit is calculated with the 



reduced density matrix: I^^ = ~Ip'Sp{pal'''). The results 
of the calculations are also presented as the color insets in 
Fig. [T2l' f) for the following parameters: the strength of dissi- 
pation a = 0.1 and the driving ampHtude /ac = 8 x 10"'^. 

Lasing in the two-qubit system 

Consider now the influence of the dissipation on the dynam- 
ics of a two-qubit system. For this the Bloch-Redfield formal- 
ism will be used. The strong dependence of the inter-level 
relaxation rates on the controlling magnetic fluxes will be 
demonstrated for the realistic system. This allows to propose 
several mechanisms for lasing in this four-level system [54]. 

For identification of the level structure and understanding 
different transition rates it is instructive to start from consid- 
ering the case of two non-interacting qubits, that is J = 0. 
In this simplified situation, the energy levels of the system of 
two qubits consist of the pair-wise summation of single-qubit 
levels, 

AE,, , 1 / . (39^ 



r(i) , 



= ±- 



,(0)2 



2 2 

In Fig. \T4[a) the energy levels are plotted as a function of the 
partial bias in the second qubit /;,, fixing the bias in the first 
qubit fa- Then the single-qubit energy levels appear as the 



horizontal lines at — 



.(0)2 



and as the parabo- 



las at E^ifh) = ±y e^°\fbf + Al For the lasing the 
hierarchy of the relaxation times is required. For this it is nat- 
ural to assume that the relaxation in the first qubit is much 
faster than in the second qubit. This allows to consider three- 
and four-level lasing schemes in Fig. Ob, c). 

As a next step, the interaction of the qubits, J 7^ 0, should 
be considered. To describe the relaxation in this system, the 
operators are converted to the basis of eigenstates of the un- 
perturbed Hamiltonian. In this representation Hq = S^^H^S 
is the diagonal matrix; the unitary matrix S consists of eigen- 
vectors of the unperturbed Hamiltonian; the excitation opera- 
tor V{t) is converted as following 



(60) 

The dissipative environment can be described as the ther- 
mostat, for which the convenient model is the bath of har- 
monic oscillators, see Fig. |9] Within the Bloch-Redfield for- 
malism, the Liouville equation for the quantum system inter- 
acting with the bath is transformed into the master equation 




FIG. 14: (Color online) Energy level structure with J = 0. (a) 
One-qubit and two-qubit energy levels as functions of the magnetic 
flux fb at fixed flux fa ■ The arrows show the fastest relaxation, which 
is assumed to relate to the qubit a. (b) and (c) Schemes for three- and 
four-level lasing at fb — fbL and = fbK- The driving magnetic 
flux pumps (P) the upper level; fast relaxation (R) creates the popu- 
lation inversion; the two operating levels can be used for lasing (L). 



for the reduced system's density matrix p(t). Then the master 
equation for the density matrix of our driven sys tem can be 
written in the energy representation as following 



) Pij ' 



Sij^PnnWjn-JijPij, (61) 



where ujij — {Ei — Ej)/h, and the relaxation rates W„ 
2ReTnmmn ^nd 



nnmm 



(62) 



are defined by the relaxation tensor Timnk, which is given 
by t he Fermi Golden rule. As it was shown in Refs. 11147 - 
I149I1 . the noise from the electromagnetic circuitry can be 
described in terms of the impedance Z{ljj) from a bath of 
LC oscillators, described by the Hamiltonian of interaction 

i/i = i ^cri"'' + crp-'^ X in terms of the collective bath co- 
ordinate X = X^fe'^fc'^fe- Here $fc stands for the magnetic 
flux in the fc-th oscillator, which is coupled with the strength 
Ck to the qubits. It follows that the relaxation tensor Timnk is 
defined by the noise correlation function 8(10) 



Timnk = ^^SM, 5H= I dte-'^'{X{t)X{0)}, 



ft2 



A 



Imnk 



(63) 



Im \ J nk 

The correlator S{uj) was calculated in Refs. fl47, 148 1 within 
the spin-boson model and it was shown that the relevant real 
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part of the relaxation tensor 

8n 



coth ■ 



huJrj 



2kBT 
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(64) 



is defined by the environmental Ohmic spectral density 
J(w) = aHu! and is cut off at some large value cjc, where 
a is a parameter that describes the strength of the dissipative 
effects. 

From the above equations the expression for the relaxation 
rates from level In) to level I to) follows 
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(65) 



InRef. (sA 

these relaxation rates were calculated as functions 
of the partial flux biases fa and /{,. It was demonstrated that 
the fastest transitions are those between the energy levels cor- 
responding to changing the state of the first qubit and leaving 
the same state of the second qubit. Such a difference in the 
relaxation rates creates a sort of artificial selection rules for 
the transitions, similar to the selection rules studied e.g. in 
Refs. [150-152]. To describe the hierarchy of the relaxation 
rates, consider them in the simplified case, ignoring the in- 
teraction between the qubits; then the single-qubit relaxation 
rates follow from Eqs. (US) and ^ 1 102, . 1531 



Re7oi = 



aA^ , AE 

— coth , 

2hAE 2T ' 



(66) 
(67) 



h AE^' 

In particular, in the vicinity of the point ~ /^* in Fig. fl4l' a). 
where AE^'''> = AE^^\ we obtain T^"^'^ /T^''^ ~ (At/Aaf. 
If Aa ^ Ab is chosen, consequently the first qubit relaxes 
much faster 

After the parametrization of the density matrix, pij 



iyij, the system's dynamics is described by the equations Il54tl 



T y]^^ + J2 ^irXrr " W^.Xu, i = 1, 2, 3; 
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0, Xoo = 1 - {Xll + X22 + 2:33); 
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When discussing Fig. [14] we pointed out that in the system 
of two coupled qubits there are two ways to realize lasing, 
making use of the three or four levels to create the population 
inversion between the operating levels. In Ref. [541 the lasing 
in the two-qubit system was demonstrated by solving numer- 
ically the Bloch-type equations ( l68b . Besides demonstrating 
the population inversion between the operating levels, an ad- 
ditional signal with the frequency matching the distance be- 
tween the operating levels was applied, to stimulate the tran- 
sition from the upper operating level to the lower one. So, the 




FIG. 15: (Color online) Two-qubit lasing and stimulated transi- 
tion. The time-dependent occupation probabilities are plotted for 
one- and two-photon driving. The driving and fast relaxation cre- 
ate the inverse population between the levels |2) and |1); then the 
stimulating signal /l cos is turned on. (54ll 



driving was considered to be, first, the monochromatic signal 
f{t) = /ac sinwt to pump the system to the upper level and 
to demonstrate the population inversion. Then another signal 
stimulating transitions between the operating laser levels is 
applied f{t) — /ac sinwt + /l sin with Scjl = E2 — Ei. 
Solving the system of equations ( |68] |. one obtains the popu- 
lation of i-th level of our two-qubit system. Pi = xu. The 
results of the calculations are presented in Fig. [15] where the 
temporal dynamics of the level populations is given for two 
situations. 

As shown in the inset schemes in Fig. [15] the fastest (dom- 
inating) relaxation transitions are |3) |2) and |1) |0). 
The system is excited by either one- or two-photon transitions, 
with fkj — E^ — Eq in Fig.[T5ja) or with 2hLj = E^ — Eq in 
Fig. fTSl b). This creates the population inversion between the 
levels 1 2) and 1 1 ) . Note that analogous competition of the driv- 
ing and relaxation can lead to the population inversion in other 
multi-level systems 1 124> .1541 . Fast relaxation, |1) — |0), 
helps creating the population inversion between the laser lev- 
els 1 2 ) and 1 1), which is the advantage of the four-level scheme 
il55[l . Then the transition |2) — > |1) is stimulated by an- 
other signal with a frequency matching the laser operating 
levels (ftwL = E2 — El). Figure [15] was calculated for 
the following realistic parameters 153]: Aa/h — 15.8 GHz, 



Ab/h = 3.5 GHz, 4°^$o//i 375 GHz, I^"'^o/h = 700 
GHz, J/h = 3.8 GHz, keT/h = 1 GHz; and also a;L/27r = 9 
GHz, /l = /ac — 5 X 10~^ with the driving frequency 
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a;/27r = 47.4 GHz for (a) and w/27r = 23.7 GHz for (b). 

For the realization of such lasing schemes, the system of 
two qubits should be put in a quantum resonator, e.g. b y cou- 
pling to a transmission line resonator, as in Ref. Il44[ Then 
the stimulated transition between the operating states, demon- 
strated in Fig. [15] will result in transmitting the energy from 
the qubits to the resonator as photons. 

CONCLUSIONS 

Here we presented the experimental and theoretical results 
of the study of driven single and coupled superconducting 
qubits. The multiphoton transitions in both charge and flux 
qubits were studied in details. Those processes are important 
for both demonstrating the fundamental quantum phenomena 
in mesoscopic systems and for developing controlling mecha- 
nisms for perspective devices. 

The system of qubits, coupled to the controlling electronics 
and measuring resonator, can be described within the semi- 
classical approach. After presenting this formalism in ap- 
plication to probing the qubit systems, we have shown some 
specific experimental results, which were accompanied by the 
calculated counterparts. The agreement between them shows 
contemporary possibility to demonstrate and describe quan- 
tum phenomena in mesoscopic systems. 
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